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Abstract. We discuss, at the mean-field level, the asymptotic shape of the reaction fronts in the general 
nA + mB C reaction-diffusion processes with initially separated reactants, thus generalizing to arbitrary 
reaction-order kinetics the work done by Galfi and Racz for the case n = m = 1. Consequences for the 
formation of Liesegang patterns are discussed. 

PACS. 82.20.Mj Nonequilibrium kinetics - 82.20. Db Statistical theories - 66.30. Ny Chemical interdiffusion 



1 Introduction 

Consider two reactants, initially separated, which are put 
in contact at time t = and start to mix one into each 
other by diffusion. A region where the reaction rate is high 
will develop at their interface. The mathematical function 
describing the variation, in space and time, of the ampli- 
tude of the reaction rate in this region is usually called a 
reaction front i?(x, t). 

Dynamical properties of reaction fronts in (purely or 
effectively) one-dimensional reaction-diffusion systems but 
that reduce to an effectively of a symmetry in the ini- 
tial of numerous studies @,|,|,|,|,|,|,||,|l^,[l|,|2|. In 
general, it is observed that these fronts obey asymptotic 
scaling, characterized by a scaling function G and scaling 
exponents a, 7 : 



R - t-'^^ 



t" 



(1) 



In the previous expression, x f (t) locates the position of the 
front (usually defined as its first moment), which generally 
obeys 



Xf{t) (X Vt 



(2) 



accounting for the diffusive origin of the front's dynamics. 

In the framework of a mean-field approximation (which 
is ours from now on), scaling hypothesis, together with 
balance considerations, can lead quite directly to the val- 
ues of the dynamical scaling exponents (see 0,^). Ac- 
cessing the structure (i.e. the shape) of the scaling func- 
tion itself requires however to go one level deeper into the 
analysis of the process in consideration. 

In addition to satisfying a purely theoretical curiosity, 
knowing better itself can provide a practical advantage. 
For example, in studying phenomena involving reaction- 
diffusion processes, like Liesegang patterns formation, it 



may be useful to dispose of an explicit, approximated an- 
alytical form for This allows, for example, to bypass 
the dynamical generation of this front in a numerical sim- 
ulation and save computation time [p^ . Finding such an 
analytical approximation requires evidently to gain suffi- 
cient information about the scaling function. 

Both tasks (derivation of the scaling exponents and of 
the scaling function) were accomplished in the pioneer- 
ing paper by Galfi and Racz where they studied the 

k 

reaction front in the A + B ^ C process with initially 
segregated As and B-s in mean-field approximation. 

In the present paper, we provide a generalization of 
their work to the case of arbitrary reaction-order kinetics, 

k 

nA+mB — > C, and explain how to calculate the associated 
C density profile in the asymptotic regime. 

An important motivation for this generalization is the 
following : in the case of Liesegang patterns, the primary 
chemical reaction leading (through several complex coars- 
ening processes) to the formation of precipitate turns out 
to be most often of the types A + 2B or 2 A + B, and not 
A + B, as usually considered for simplicity in theoretical 
models (l3[|l|,|l|. 

2 Scaling analysis. 

2.1 Definition and notations. 

The case of mean-field, general reaction-order kinetics in 
the initially segregated reactants case has already been ad- 
dressed in [^. Using scaling analysis, the authors showed 
that the exponents controlling the asymptotic behaviour 
of the reaction front are given, in terms of the reaction- 
order constants, by : 



a(n, m) 



n + m — 1 



2{n 
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where a, 7 are the same as in (y), and (n, m) are integers The equations then read 
both > 1. 

It is important to note the foUowing properties of these 
exponents : 



1. a{n, m) = a{n + m) , j(n,m) = j(n + m). 

2. a{n,m) < 1/2 \/ {n,m), a increases monotonicany 
from 1/6 (n = m = 1) to 1/2 (n + m — > 00). 

3. a{n,m) +j{n,'m) — 1/2. 

We can start from the above resuhs to formulate a gen- 
eral derivation that will lead us to the family of ordinary 
differential equations defining the asymptotic shape of the 
reaction fronts Q- To this goal, let us consider the 

following one dimensional initial- value problem, describ- 
ing a reaction-diffusion process between initially separated 
A and B particles in the mean-field approximation : 

drAiX, T) = DAd\A{X, T) - kn{A''B"'){X, T) (4a) 
dxBiX, T) = DBd\B{X, T) - km{A'^B"'){X, T) (4b) 
dTC{X,T) 

with 



= fcA"(X,T)B'"(X,T), 


(4c) 


A{X,T = Q) = aoe{-X) , 


(5a) 


B{X,T = Q) = bQe{X) , 


(5b) 


C{X,T = 0) = 0. 


(5c) 



In the above equations, 

— 6 denotes the Heaviside step function [9{X < 0) — 
0, e{X>Q)^l]. 

— A, B and C are concentrations with dimensions 
[A,B,C]=X-\ 

— Da and Db are diffusion coefficients {[Da,Db] = 
X^T-i). 

— fc is the reaction rate {[k] = [X'^+"^-^T-^]). 

In the following, we will only consider the case of equal 
diffusion coefficients, Da = Db = D, since the method we 
are going to use requires this strong condition to be satis- 
fied. The asymmetric case Da 7^ Db reveals to be several 
orders of magnitude higher in difficulty. Interesting results 
have been obtained in the case n = m = 1, in connection 
to the front's dynamics but a derivation of the shape 
of the scaling functions for arbitrary Da/ Db seems to be 
still out of reach for the moment. 



with 



dta{x, t) 


= d^a{x, t) 


-na"(a;,i)6™(x,t) 


dtb{x, t) 


^dlb{x,t) 


~TOa"(x,t)fe'"(a;,i) 


dtc{x,t) 


= a"{x,t)b'''{x,t) , 




a(x, t = 0) 


= 9{-x) , 




b{x,t = 0) 


= Nx) , 




ao 




c{x,t^ 0) 


= . 



2.2 Solution for a — {n/m)b. 
We define 

u(x, t) = (a 5 J (x, t). 

\ ml 

This function obeys the diffusion equation : 

dtu{x,t) ~ d'^u{x,t) , 

u{x<0,t = 0) = l , 

n bn n 

u{x > 0,i = 0) = = q 

m ao m 



whose solution reads 

1 n / n ^ „ 



(7a) 
(7b) 
(7c) 



(8a) 
(8b) 
(8c) 



(9) 



(10a) 
(10b) 

(10c) 



(11) 



In the above equation, erf denotes the error function, 
erf (a;) = (2/0?) exp{—w'^)dw. 

Let Xf{t) be such that u{xf{t),t) — 0. One can check 
that 



with Df = Df{q) given by 




m. ^ 



(12) 



(13) 



The first step in our calculation is to render a-dimen- 
sional the problem we are dealing with. This can be done 
through the following change of variables : 



V D 



X 



a,b,c = A/ao,B/bQ,G/co 



(6a) 

(6b) 
(6c) 



^ This derivation follows closely the steps and notation in [Q] , 
and the reader should refer to it for further details and justifi- 
cations. 



2.3 Equation for a in the reaction zone. 

We write now b = ^{a~ u) and plug it into ( [Za] ) , getting 
thereby an equation for a involving only a and the known 
function u : 



dta{x, t) = d^a{x, t) — n y — j 



a [a — u) 



ix,t) 



(14) 



We are interested in the solution of this equation in the 
reactive region |a; — « ia{n,m) ^ ^^^^ latter is believed 
to widen with a time exponent a{n, m) < 1/2 , this allows 
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us to expand u around xj to the lowest-order in x/^/t, 
since the neglected terms will vanish as t — > oo : 



Vt 



with K given by 



l + -q 



(16) 



The boundary conditions that the solution to (Q) must 
satisfy in the reactive region are : 



a{x — !■ —GO, t) = —K 
a{x — > -|-oo, i) = 



(17a) 
(17b) 



2.4 Scaling hypothesis. 



We shall now assume that asymptotically (i.e. when t — > 
oo), the solution to (llj) adopts the following scaling form : 



a(a;,i)«t-^("'™)G(„,™) 



(18) 



where {G(^n.m)}n.m>i are a family of scaling functions re- 
maining to be characterized. The scaling exponents are 
given by (|). 



2.5 Differential equation for G'(„ m)- 

Let's define first the reaction zone coordinate z 

X-Xf 



(19) 



Inside the reaction zone, u and h = — ia — u) write 

u[z) = -^^("^'"'-i/^z , (20) 



6(z) = -(i-''("'™)[G(„,,„)(z)-f Xz]) . 

Using (|l8|), eq. (^ becomes : 
i'"^"'"^-'[-7(^,HG(„,™) - a(n,™)za,G(„,™)] 



(21) 



We now take the asymptotic limit inside the reaction 
zone, i.e. we let < — > oo, keeping z fixed. The two terms on 
the left-hand side vanish (remember that a(„^m) < 1/2 !) 
and we remain with the following ordinary, non-linear 



second-order differential equation for the scaling functions 
G 



(n,m) 



(15) G'('„,™)W = "(-) Gf„_^)W[G(„.™)(^) 



Kz] 



(23) 



The boundary conditions (17a-17b|) imply the follow- 



ing asymptotics for G 



(n,m) 



G{n,-m){z) ^ -Kz, Z ^ -CO 
G(„,m)(z)^0, 2 ^ OO . 



(24a) 
(24b) 



We are now left with a boundary- value problem (^ 24 ) 
that can be solved numerically. 



2.6 Solving the equation for G(„ m). 

We can make the problem K-independent by rescaling G 
and z : 



G = i^^("'™)G , 
z = K^'^'^'^'h , 



(25a) 
(25b) 



and by using a suitable choice for fi and i^. Inserting these 
scaled forms into ( p3| ) and imposing that K drops out leads 
to 



/i(n, m) 



n + m + I 



= ^i{n + m) 



-{n + m- 1) 

v(n,m) = = i/in + m) 

^ ' ' n+m+1 ^ ' 

The problem we are left to treat is now : 
GlJ~z) = n (-) Gl„,Cz) [g„MS) + S 

Gn,m{z) -Z Z — > OO , 

Gn,m{z) ^0 Z ^ OO . 



(26a) 
(26b) 



(27a) 

(27b) 
(27c) 



The reader should keep in mind, from now on, that the 
"tilde" sign stands for quantities expressed in terms of the 
rescaled, if— independent version of the scaling function 
G(n.m) and reactive coordinate z. 



2.7 The dimensionless reaction front. 

By definition : 



(22) i?(„,™)(x,t) = a"(x,t)6"(x,t) 



t-0MF^^^.^){z) , 



Kz] 



(28) 



the last inequality defining both the reaction rate ampli- 
tude exponent /3 and the asymptotic reaction front scaling 



function Fi 



{n,m) • 
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S. 0.2 



(n,m)=(1,1 

(n,m)=(1,2) 
- (n,m)=(2,1) 
(n,m)=(2,2) 




Fig. 1. Numerical solutions to ( ^t\) for (1 < 



<2). 



Figure |i| shows the resuh of the numerical computa- 
tion of the fronts F^^.m) = F{G, z) for the cases 2 < 
n + m < A (only such low values of n -|- m are relevant 
in connection to experiments). The reader should not be 
surprised by the asymmetry between the {n,m) = (1,2) 
and (n, m) = (2, 1) cases : it is due to the passage from 
the variables {G(„.„i),z} to {G(„^m),z}, since both quan- 
tities are rescaled by a K— dependent factor which is not 
(n, m)-symmetric ! 

These fronts are all essentially localized around z — 0, 
as Fig. |] suggests. In fact, one can check that in the sim- 
plest (n,m) = (1,1) case, the dominant contribution to 
both tails at z ^ ±oo is proportional to exp(— z^/^). In 
the cases n = l,m > 1, the decay can be shown to be 
dominated by exp(— z^+™/^) at z ^ oo but only alge- 
braically for z ^ -cx) : i?(z -> ~oo) w ^-(2m+i)/(m-i)_ 
Finally, when n > 1, m > 1, both tails are algebraic, with 

R{z oo) « ^-(2™-Kn)/("-l)^ 

It is also worth noting that in the asymmetric cases 
n ^ m, Xf does not coincide with the maximum of the 
front. 



3 The C concentration profile. 
3.1 Derivation of the asymptotic profile. 

We are interested now in estimating the (possibly x-depen- 



dent) density 



x) of C particles left behind by the 



fronts R(n,m) which travel diffusively through the system. 
This quantity is, for example, of great importance in the 
theories of Liesegang pattern formation |jlJ,Q. 

In dimensionless units, Cq"'™'' is formally given by 



(n.m) , 

Co ■ ' 



(x) = / R(ri,rn){x,t)dt 



(29) 



(n.rn) 
Co 

of 



[x) turns out however to be possible by making use 



a) the precious algebraic relation a + ^ = 1/2 between 
the scaling exponents, and 

b) the particular structure of the solutions to ( ^7|) . 

Let's consider first a narrow slice SF(^n „i){zo,S'w) of 
the scaling function F(^n,m)i centered on z = zq, of width 
Sw <C 1. In the spirit of the Riemann integral, we can 
approximate the amplitude of 5F inside [zq — 5w/2, zq + 
Sw/2] by its value i^(„,m)(-2o) at the center. 

We can estimate the contribution SC^'^'"^\zo,x)/6x of 
this slice to the asymptotic local C density inside [x, x+Sx] 
as follows. The fraction of the front we are considering 
will reach a; at a certain time t(zQ,x). The quantity of 
C particles deposited in the interval [a;, a; -I- Sx] will be 
proportional to the amplitude, the width and inversely 
proportional to the speed of the slice at t = t{zo,x) : 



5C("^"')(zo,a;) w 

t(zo,a;)-^("'"')F(„,^)(^o)t(zo,x)°("'"')^«; 
^Df/{2t{zo,x)) 



6x 



^2lDfFi^n^„,) {zo)6w 5x 



(30) 



In other words, the contribution of the slice to the 
density at x is proportional to its "mass" {zq)5'w but 

independent of t(zo, x), and hence of a;. This indicates that 
the asymptotic density profile is fiat. By superposition, our 
argument leads immediately to the result we are looking 
for : 



(a:) 



(n.m) 



const. 
2/Df I F^„^^){z)dz 



(31) 



Now our real fortune is that we are able to evaluate 
analytically J F(n.m)- Using (|2^ ) and (p8|), We have 



1 



5^G(„^„j)(z -oo) 



(32) 



It is intuitively clear, from purely physical considera- 
tions, that the solution will converge to its values at ±oo 
in such a way that 

lim G"(„.,„) (z) = -X , lim G(„ ,„)(z) = (33) 



[see (p7q)-(27c)]. So we finally have 



K 

F(n,-m){z)dz = — 



(34) 



and we end up with the result 



Due to the several timescales dependence of R(n,m)j 
this integral is difficult to handle. The estimation of 



(n,m) 
^0 



(35) 
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Going back to the dimensional variables A, B, C, X 
and T, one can check that (351) writes 



(n,m) 
-0 




SO we obtain from (|37|) 



(36) 




e^p{-Df/2D) 



1 D 
2D] 



(39) 



l + -q 

which, together with the dimensional expression for K, 
1 + -q 

in. ^ 







V___(n,m)=(2,1) 


0.01 




(n,m)=(1,1) 


f 






cf 

o 
o 




(n,m)=(1,2) 


0.005 
n 







K 



■eyi^{~Df/2D) 



(40) 



gives finally 



(n,m) 
•^0 



m 



ID ,D 

+ — f 

2Df \}Df' 



The physical meaning of this result is clear: if Da = 
Db = D and g ^ 1, then Df ^ D and the B particles 
appear as nearly immobile for the invading As. As m B-s 
are required to produce one C, the density equals b^/m to 
a very good approximation. 

However, in the typical conditions of a Liesegang ex- 
periment ( where Da — Db and 10^^ < <? < 5 • 10"^ 



-500 



500 1000 
X 



1500 2000 typically), the first order correction in Z3/Z)y to 



(n.m) v 

3A ■ ' lies 



Fig. 2. Numerical solutions to (^j for (2 < ti + m < 3) (plain 
curves). The dashed lines indicate the asymptotic values as 
calculated from (|l^), (|l^) and (|3^. The following values of 
the "free" parameters were chosen for all three cases -.0 = 1, 
ao = 10~^6o = 1, fc = 0.1. Units for X and k vary with (n, m) 
according to (E3), while Cq"'™' is measured in units of l/ao- 



Fig. |2| shows the Cq '™ profiles obtained by numeri- 
cal integration of the reaction-diffusion equation s (0 ), to- 
gether with the asymptotic values predicted by (p6]). 



3.2 Low q expansion. 



in the range (0.1 — 0.2){bo/m), and should therefore, in 
principle, not be neglected, as can be seen of Fig. ^j. 



4 Summary. 

We have derived the family of ordinary differential equa- 
tions defining the asymptotic shape of the reaction fronts 

k 

in the nA + mB — > C reaction-diffusion process with ini- 
tially separated reactants [eqs. (p3|),(24)]. The four lowest- 
order cases in n+m have been solved numerically (Fig. Q). 
We have also shown, and confirmed by numerical simula- 
tions, that the density Cq"'™^ of C particles deposited in 
the system by these traveling fronts is asymptotically con- 
stant (Fig. 0), and we have made explicit the dependence 
of this density on the reaction orders n, m, as well as on 
the material parameters D, aq, fc and bo entering the prob- 
lem [eq. (|36|)]. 



In the context of Liesegang patterns- forming experiments, 
which are generally based on the penetration of a highly 
concentrated solution into a dissolved one, it is useful to 
dispose of an expansion of this result for low q = bo/ao 
values. We recall that D/ is given by 



n^ = l-2-q - 1 . (37) 



The large x asymptotics of Erf(a:) is given by ||n 

9 

1 




Erf(a;) 



1 - 



1 - 



2a;2 



(38) 



5 Conclusion. 

To conclude this study, we would like to comment on the 
interesting phenomenon shown by Fig. |2| : when, for the 
two TO = 1 cases, no significant quantity of reaction prod- 
uct is created in the majority species subspace, the case 
where m = 2 exhibits, on contrary, an important deposit 
of C on the left hand side, up to far beyond the initial 
location of the interface. This fact must be evidently re- 
lated to the details of the short-time dynamics of the re- 
action front. Some studies have already been carried on 
the early-time regime subject for the n = m = 1 case 
[Tsl in the past. They unveiled the existence of a surpris- 
ingly complex behaviour, including successive power-law 
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regimes for the early front's dynamics, and even the pos- 
sibiUty of a change in its direction of motion. Such non- 
trivial behaviour has also been observed numerically in 
the higher-order kinetics cases we have addressed in the 
present paper, and a detailed study of the dependence on 
(n, m) of the short-time dynamics should be the object of 
a forthcoming paper. 
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